Quantum theory of resonantly enhanced four-wave mixing: mean-field and exact 

numerical solutions 
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We present a full quantum analysis of resonant forward four-wave mixing based on electromagneti- 
cally induced transparency (EIT) . In particular, we study the regime of efficient nonlinear conversion 
with low-intensity fields that has been predicted from a semiclassical analysis. We derive an effec- 
tive nonlinear interaction Hamiltonian in the adiabatic limit. In contrast to conventional nonlinear 
optics this Hamiltonian does not have a power expansion in the fields and the conversion length 
increases with the input power. We analyze the stationary wave-mixing process in the forward 
scattering configuration using an exact numerical analysis for up to 10 3 input photons and compare 
the results with a mean-field approach. Due to quantum effects, complete conversion from the two 
pump fields into the signal and idler modes is achieved only asymptotically for large coherent pump 
intensities or for pump fields in few-photon Fock states. The signal and idler fields are perfectly 
quantum correlated which has potential applications in quantum communication schemes. We also 
discuss the implementation of a single-photon phase gate for continuous quantum computation. 



I. INTRODUCTION 

The cancellation of linear absorption and refraction 
in resonant atomic systems by means of electromagnet- 
ically induced transparency (EIT) [Q, [| has led over 
the past 10 years to fascinating developments in nonlin- 
ear optics [Q, |). For example, coherently driven, res- 
onant atomic vapors under conditions of EIT allow for 
complete frequency conversion in distances so short that 
phase matching requirements become irrelevant ||. It 
has been predicted that resonantly enhanced nonlinear 
interactions of light and atoms based on EIT will lead 
to a new regime of efficient nonlinear optics on the level 
of few photons jj], ||, [|. Besides being of interest in 
its own right, such a regime would be very important for 
applications in quantum communication and information 
processing. It is clear that quantum effects will play an 
essential rule in this regime and the quantum dynamics 
may substantially deviate from semiclassical predictions. 
With the exception of a few exactly solvable problems, 
for example the resonantly enhanced Kerr effect || |To|j . 
quantum treatments of EIT-based nonlinear optics have 
so far been restricted to small-fluctuation approxima- 
tions. In this paper we present a full quantum analysis 
of a particular EIT-based nonlinear system, namely reso- 
nantly enhanced four-wave mixing in a doublc-A system 
with co-propagating pump modes [jll| [l2| [f| [lj], |l5) . 

Within a semiclassical analysis it has been shown that 
resonantly enhanced four-wave mixing can lead to com- 
plete conversion of the pump-field energy into the signal 
and idler modes even for very weak pump fields. If the 
atomic degrees of freedom can be eliminated adiabati- 
cally and losses can be ignored, the semiclassical nonlin- 
ear problem is exactly integrable |l6|, For counter- 
propagating pump modes a phase transition to mirrorless 
oscillations has been predicted |l8|, [ll], and experi- 
mentally verified f2l| . A linear fluctuation analysis has 
shown that close to the threshold of oscillation an almost 
perfect suppression of quantum noise of one quadrature 



amplitude of a combination mode of the generated fields 
occurs [^2], In addition, sufficiently above thresh- 

old, light fields with beat-frequencies tightly locked to 
the atomic Raman-transition and extremely low relative 
bandwidth are generated p6[ . 

Assuming conditions of adiabatic following and con- 
sidering the limit of an infinitely long lived ground state 
coherence we here derive a classical effective nonlinear 
Hamiltonian which only contains field degrees of freedom. 
In contrast to conventional four- wave mixing |2^, ^3|, ^4) , 
this Hamiltonian is a ratio of polynomial expressions 
and has no power expansion in the fields. As a conse- 
quence the conversion length increases rather than de- 
creases with growing input power. The evolution corre- 
sponding to this classical Hamiltonian can be mapped to 
a single nonlinear pendulum and can be solved exactly. 
However, the initial state of the pendulum correspond- 
ing to vacuum in signal and idler modes is an unstable 
equilibrium point. Thus the initial evolution is entirely 
governed by quantum fluctuations. Replacing the classi- 
cal field variables in the effective Hamiltonian by oper- 
ators in normal ordered expressions, we obtain a quan- 
tized Hamiltonian. Due to its nonpolynomial character it 
is not possible to apply phase space techniques to study 
the quantum evolution of the fields. Instead, making use 
of constants of motion, the stationary four-mode interac- 
tion is reduced to a single-mode problem, which can be 
solved numerically for up to 10 3 input photons. 

We find the quantum dynamics to be significantly dif- 
ferent from the semiclassical prediction. In particular, 
complete conversion is achieved only for input fields in a 
few-photon Fock state or asymptotically for a very large 
coherent pump. The main features of the dynamics, 
such as conversion efficiency and the dominant oscilla- 
tion frequency are reproduced by an appropriate mean- 
field theory which takes into account anomalous corre- 
lations. Finally the quantum statistics and correlations 
of the fields are analyzed and potential applications for 
quantum communication and information processing dis- 
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cussed. For example, we show that the resonant four- 
wave mixing process is an excellent source of quantum 
correlated photon pairs and can be used as a single- 
photon phase gate for continuous quantum computation. 



II. ELIMINATION OF AC-STARK INDUCED 
NONLINEAR REFRACTION AND DERIVATION 
OF AN EFFECTIVE HAMILTONIAN 

The standard resonant four-wave mixing scheme in a 
double-A system is shown in Figure [l]. 




FIG. 1: Atoms in a double-A configuration interacting with 
two driving fields and two generated fields (£1,2). 



The interacting beams consist of two pump fields with 
frequencies Vdi and 1^2 and Rabi-frequencies f2i and f2 2 , 
and two generated fields (signal and idler) described by 
the complex Rabi-frequencies E\ and E2, with carrier 
frequencies v\ = v^i + and v-i = v&% — loq, where 
luq = L02 — Wi is the ground-state frequency splitting. It 
is assumed here that the pump and generated fields are 
pairwise in two-photon resonance and thus in four-photon 
resonance. The latter is a consequence of energy conser- 
vation. The two photon resonance is a consequence of 
phase matching and semiclassical treatments show that 
signal and idler fields are generated precisely at those 
frequencies. The fields interact via the long-lived co- 
herence on the dipole-forbidden transition between the 
met ast able ground states |1) and |2). 

The problem with this model is that associated with 
the finite detuning A, which is necessary to minimize 
linear absorption, are ac-Stark induced nonlinear phase 
shifts. These phase shifts reduce conversion efficiency 
from the pump to the generated modes, and at the same 
time increase the distance required for conversion to take 
place [^?J. As was shown in |27) these problems can be 
overcome by modifying the system slightly. Instead of 
the original four-level scheme a five-level set-up depicted 
in Figure ^| is used. This symmetric configuration cancels 
ac-Stark shifts. In order to maintain the nonlinear inter- 
action, which is also an odd function of the detuning, it is 
however necessary to choose atomic states such that the 
coupling constant for one of the four transitions |1) — > |3), 
|1> -> |4), |2) -> |3), and |2) -> |4) has a different sign 



to the other three. This can easily be accomplished by 
using different hyperfine levels |27| . 




FIG. 2: Modified double-A system. Addition of a fifth level 
allows the cancellation of destructive phase-shifts. 

Choosing the driving field f2i to be in resonance with 
the 1 2) — > 1 5) transition while the second driving field SI2 
has a detuning ±A with |A| 3> \^2\ ensures that the lin- 
ear losses due to single-photon absorption are minimized. 

The interaction of the fields with an ensemble of five- 
level atoms is described by Maxwell's equations for the 
fields and a set of density matrix equations for the atoms 
which determine the atomic polarizations. Under condi- 
tions of adiabatic following and for the case of classical 
fields the latter are ususally solved in steady-state and 
feed back into the first, yielding nonlinear field equations 
of motion. This procedure is rather involved, particularly 
if several atomic levels need to be taken into account, and 
is completely inadequate for quantized electromagnetic 
fields. Consequently in this paper we adopt a different 
procedure. We first derive an effective classical interaction 
hamiltonian in the adiabatic limit and then quantize the 
effective theory by replacing the field amplitudes by prop- 
erly ordered operator expressions. 

In order to rigorously calculate the medium response 
to classical fields, one would have to solve the atomic 
density matrix equations to all orders in all fields tak- 
ing into account all relaxation mechanisms. Instead we 
use a simplified open-system model which allows to de- 
rive rather compact expressions for the atomic suscep- 
tibilities. In this model the interaction of an individual 
atom with four classical modes is described by a complex 
Hamiltonian, which, in a rotating-wave approximation 
corresponding to slowly-varying amplitudes of the basis 
(|1) |2) |3) |4) |5)) T , can be written as 
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At the input the signal and idler modes E± and E2 are as- 
sumed to have zero amplitude and all atoms are in state 
|1). Taking into account that optical pumping out of |1) 
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is negligible as A 3> IO2I, this is a stable configuration 
corresponding to an approximate adiabatic eigenstate (to 
lowest order in A -1 ) of H . We now assume that the in- 
teraction takes place over sufficiently long time scales, 
ensuring that the atoms always stay in this approximate 
adiabatic eigenstate. Under these conditions the Hamil- 
tonian can be replaced by the corresponding eigenvalue 
Ao . Solving the 4th order characteristic equations for the 
eigenvalues of ([[]) and expanding the solutions in a power 
series of (51/ A), with SI being a characteristic value of the 
Rabi-frequencies, one finds to lowest order 



stants of motion: 



An = -r 



fid 



\Ei\ 



(2) 



One recognizes two interesting and unusual features. 
First the eigenvalue and correspondingly the medium po- 
larization cannot be expressed as a polynomial in the field 
amplitudes. Thus the resonant interaction corresponds 
to an all-order nonlinear process. Secondly despite the 
resonant interaction An has no imaginary component and 
there are hence no linear losses. This is due to quantum 
interference associated with EIT. It allows for efficient 
nonlinear interactions close to resonance without suffer- 
ing from linear absorption. 

To quantize the interaction problem we replace the 
complex amplitudes of the Rabi-frequencies in (^|) by pos- 
itive and negative frequency components of correspond- 
ing operators, choosing normal order in the numerator 
and the denominator, multiply by the density of atoms 
and integrate over the interaction volume. We thus ar- 
rive at the effective interaction Hamiltonian between the 
quantized electromagnetic fields 



Hi* 



hNA 



dz 



E[Ei 



(3) 



the number density of atoms, A 
cross section of the beams, and 



where N is 

the effective 

E x {z,t) = dsi E fc V"i/ 2h£ oVa lk {t)e-^-^ z / c 
is the slowly varying positive frequency opera- 
tor of the signal Rabi-frequency. Correspondingly 
filM) = d 52 J2 k V"di/2he Vb lk (t)e-^-^^ c de- 
notes the slowly varying positive frequency operator of 
the first pump Rabi-frequency. The operators a M and 
bfj, obey harmonic oscillator commutation relations. V 
is the quantization volume, which shall be identified 
with the interaction volume and dij is the dipole matrix 
element of the \i) — ► \j) transition. It should be noted 
that numerator and denominator commute and thus 
there is no ambiguity with respect to the ordering of the 
two terms. 

One can easily verify that there are four independent 
quantities that commute with H lrA and are therefore con- 



9\Q,i + E[Ei = 



El En = 



h\Ci\EiE 2 + Ci^eIeI 



constant, 
constant, 
constant, 
constant. 



(4) 
(5) 
(6) 
(7) 



These are the quantum analogs of the classical Manley- 
Rowe relations which express energy conservation in the 
system plus an equation expressing the conservation of 
the relative phase between the fields p3fl . The existence 
of these constants of motion will considerably simplify 
the analysis. 

It should be noted that the existence of these constants 
of motion are not artifacts caused by using an effective 
rather than the full Hamiltonian. We have also pursued 
a more rigorous derivation, which involved writing the 
Heisenberg equations of motion for the atomic and field 
subsystems seperately and including decay terms. In the 
limit 7 < A and negligible decay from the two ground 
states this approach also yields (Q) - (||). Unfortunately 
the equations of motion obtained in this way have a num- 
ber of unpleasant features and are not ameanable to an- 
alytic solution, which is why we have chosen to use the 
effective Hamiltonian (||) as the starting point of our dis- 
cussion. 



III. CLASSICAL SOLUTIONS FOR FORWARD 
FOUR- WAVE MIXING 

To obtain classical solutions we use (||), and note that 
the polarization P of the medium for the probe transi- 
tions can be expressed as a partial derivative of the av- 
erage single-atom interaction energy with respect to the 
electric field or, respectively, the Rabi-frequencies Ei 



Pi = -'- 



Ndi /dHi 



dE* 



i=l,2 (8) 



A similar expression holds for the drive field polariza- 
tions. Here (...) denotes quantum-mechanical averaging, 
di the dipole matrix elements of the corresponding tran- 
sition and N is the atom density. Hence one can directly 
obtain the stationary field equations in slowly-varying 
amplitude and phase approximation: 



dEj 
dz 



Vi / dH m 



dE* 



(9) 



where r\i — Ndfuji/ (2hceo). This approach allows us to 
obtain the polarizations and equations of motion for the 
fields without calculating the density matrix . 

Using (^|), going to a comoving frame via the transfor- 
mation (z, f) — > (C = z,t = t — z/c), assuming that all 
dipole moments (averaged over orientations of the atoms) 
are approximately the same and introducing the common 
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coupling coefficient n = N d 2 u>i/(2hceo) — 3N\ 2 j/8ir, 
gives the following equations of motion for the fields: 



~1K ~ j I XU ) 



d_ 

W E2 = < 
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d( A(|ni| 2 + |£iP) 2 

^ = 



' K AprFTW)' 



A(\Q^ + \E^) 



(11) 
(12) 

(13) 



Without sacrificing the underlying physics we can assume 
that the intial intensities of the pump fields and the intial 
intensities of the generated fields are equal: 



l^il 2 = |fi 2 



|£il 2 H£ 2 | 2 . 



(14) 



As can easily be seen, the constants of motion imply that 
asymmetric initial conditions lead to exactly the same 
dynamics with the intensities merely shifted up or down 
by a constant. 

Introducing a normalized intensity y(t) with the iden- 
tification ( — > t 



m 



l^iWI 5 



ittii 2 + i£ir 

k Im(Oi0 2 -E'i^2" 



I -Si I 



(15) 
(16) 



and making use of the constants of motion (Q) - (||) to 
reduce the problem to one variable, we obtain the differ- 
ential equation 



(A/k) 2 d 2 y 
2 dt 2 



1) 



(17) 



In this form it is clear that the normalized intensity of 
the drive field can be identified with a particle of mass 
m = A/k moving in the potential V(y) = —Ay 2 (y — l) 2 , 
depicted in Figure ||. Conservation of energy requires 
that the particle be trapped between < y < 1. If 
-Ei(O) = £2(0) = 0, i.e. there is no seeding of the gen- 
erated fields, then we have the case where y(0) — 1. As 
is clear from the potential diagram, this corresponds to 
a point of unstable equilibrium, and indicates that no 
dynamics can take place. 

Choosing to seed the generated fields with a specific 
amplitude and phase corresponds to moving the particle 
off the critical point in the potential diagram. The drive 
field intensity will then oscillate, transferring energy back 
and forth between the drive and generated fields accord- 
ing to the constants of motion. The dynamics is however 
sensitive to both the strength of the seed fields and the 
initial relative phase between all four fields. The exact 
analytical solution shows that the oscillation period of 
the energy transfer as well as the maximum conversion 
efficiency depend on the seed-field intensity and phase 

MM- 




All energy in 
pump fields 



FIG. 3: Nonlinear pendulum potential experienced by a "par- 
ticle" corresponding to the exchange energy between the 
pump and generated fields. k/A taken to be unity. 



IV. QUANTUM THEORY OF STATIONARY 
FORWARD FOUR- WAVE MIXING 



A. Exact numerical solution to the stationary 
quantum problem 

While forward four-wave mixing is relatively well un- 
derstood in the semiclassical case, very little work has 
been done on the fully quantum case, that is, where 
both the atomic subsystem and the interacting electro- 
magnetic fields are quantized. There are at least three 
reasons for considering such a formulation. 

First, from ( |Io| ) - ( |l3|) we see that the nonlinear inter- 
action of the system actually increases as the strength of 
the pump fields is reduced. However, when the strength 
is reduced to a level corresponding to only a few photons, 
the semiclassical approximations must fail. Nonetheless, 
the extreme nonlinearities experienced in conditions of 
very weak fields opens the intriguing possibility of non- 
linear effects at the few-photon level. 

Secondly, a semiclassical analysis is not able to cope 
with the physically realistic situation where initially only 
two fields are present. Vacuum seeding through quan- 
tum fluctuations cannot be described in the semiclassi- 
cal framework, and one must always introduce at least a 
third field into the analysis so that conversion can take 
place. 

Thirdly, the behaviour of the system, at least in a 
semiclassical analysis, is extremely sensitive to variations 
in intial phase and amplitude of the pump and seed 
fields. When the fields are treated quantum mechani- 
cally, these parameters will often be indeterminate, par- 
ticularly when starting from a vacuum field. 

While an analytic solution of the fully quantum case 
appears intractable, it is nevertheless possible to obtain 
exact solutions numerically under stationary conditions. 
The stationary spatial evolution of the fields can be de- 
scribed by considering the time evolution of four har- 
monic oscillators 0,1,0,2 and 61,62 corresponding to the 
generated and pump fields respectively, interacting via 
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the nonlinear Hamiltonian 



hue 



b\b\aia,2 + a\a\bib 



b\bi + a\d\ 



(18) 



bi and 



with the following identification — ► a;, f2j 
k = 3iVA 2 7/87r. 

The technique we use is to describe the possible states 
of the system in a number state basis of the fields, use 
the effective Hamiltonian dl8| ) and solve the resulting 
Schrodinger equation numerically. We will use the no- 
tation njij riEx nE 2 ) for our basis vectors, where 
ROi is the number of photons in the f^i mode, uq 2 is the 
number of photons in the 2 mode and so on. 

The difficulty with this technique is that even though 
we have eliminated the atomic degrees of freedom by us- 
ing an effective Hamiltonian, because we are considering 
a four- wave mixing process the size of the Hilbert space 
will scale as n 4 , where n is a characteristic number of 
photons in each of the fields. Consequently, if we wished 
to consider, say, 100 photons in each of the beams, our 
Hilbert space would be 10 8 dimensional, and the problem 
would involve diagonalizing matrices with 10 16 elements. 
Thus, as it stands, this approach is not computationally 
feasible. 

The scaling problem can be avoided by using the con- 
stants of motion (^) - (||). Taken together, these relations 
allow us to reduce a problem with four degrees of freedom 
to just one, which is essentially the energy transfer from 
one field to another. For example, (g) states that when 
a photon is annihilated in the f^i mode, another must 
be created in the E\ mode. Reduction of the problem to 
one with a single degree of freedom allows us to choose 
the basis 



^nm = \n\ — n n-2 — n n.3 + n + n) 



(19) 



where m, 112, and 714 denote the number of photons 
initially in each of the four modes, and n is the single 
degree of freedom denoting how many photons have been 
transferred out of the pump mode. 



B. Intensity evolution and quantum limitation of 
conversion efficiency 

We first consider the evolution of the fields when the 
initial states of the pump fields are given by number 
states. This situation is computationally easy, and we 
have calculated results for initial states consisting of up 
to several thousand photons. A number state consisting 
of several thousand photons, however, is not particularly 
experimentally realistic, and thus here we consider onl 
low photon numbers. The results are shown in Figure 

We note that the solutions are oscillatory, in quali- 
tative analogy to the semiclassical predictions. In the 
fully quantum case, however, we see that only when the 
initial pump fields contain one or two photons are the 
oscillations comprised of just one frequency. If the initial 




5 10 15 20 25 30 

Interaction distance (A/k) 

FIG. 4: Evolution of average pump photon number from ini- 
tial number states. From top to bottom the initial photon 
number in the pump modes are 1, 2, 3, 5 and 15 photons. 



fields contain three or more photons several oscillation 
frequencies are present, leading to a reduced conversion 
efficiency. Provided long enough interaction distances are 
considered, however, the different oscillation frequencies 
come back into phase and reinforce each other, leading 
to the conclusion that complete conversion can always be 
obtained at some point. 

Another point to note is that as the intensity of the 
initial pump fields is increased, the oscillation period in- 
creases. This is in agreement with the form of the Hamil- 
tonians (^), (|l8|). The denominator is a constant of mo- 
tion and is related to the intensity of the pump field 
Thus a more intense pump field gives a smaller inter- 
action energy and consequently a longer conversion dis- 
tance. 

Next we consider the case where the pump fields are 
initially described by coherent states. This situation is 
considerably more computationally intensive, but calcu- 
lations with an average photon number of up to 1000 in 
each pump mode have been carried out. The results are 
shown in Figure |^. There are several overall qualitative 
features. 

We see oscillations on a short distance scale damping 
out over a longer distance scale to a conversion efficiency 
of approximately one third. The damping at longer dis- 
tances can be explained by considering a coherent state 
to be superposition of number states. Each number state 
has a different oscillation period, and consequently inter- 
fere with each other and get out of phase. As would 
be expected, on still longer distance scales fractional re- 
vivals are seen. As the input power is increased, the con- 
version distance (distance to the first mininum in pump 
field intensity) increases logarithmically while the con- 
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version efficiency asymptotically approaches unity. The 
scaling of the conversion distance with input power in the 
resonantly enhanced four-wave mixing scheme is exactly 
opposite to the case of ordinary off-resonant four-wave 
mixing. There the conversion distance decreases with in- 
creasing input power. 



C. Quantum correlations and fluctuations of the 
generated fields 

We now discuss the quantum fluctuations and corre- 
lations of the generated fields. It is immediately evi- 
dent from the interaction Hamiltonian that when the two 
generated fields start in the vacuum state, they will at 
all times be perfectly photon-number correlated. Only 
states with equal photon number in both modes can be 
generated. Consequently the intensity difference between 
the two fields is perfectly squeezed: 



Al 



a\a\ 



= 0. 



(20) 



To characterize the statistics of photon pairs in the two 
generated modes we have calculated the Q-parameter for 
different input states and intensities. As can be seen from 
Figure H the pair statistics remain sub-Possonian (Q < 0) 
for a Fock-state input with n < 2. For a Fock-state 
input with n > 2 the pair statistics have a sub-Poissonian 
character only for very small interaction distances and 
around the revival of the input intensity (see Fig. ||). 

For an initial coherent state the pair statistics remain 
super-Poissonian at all times. This is illustrated in Figure 



FIG. 6: Mandel Q parameter for initial Fock states with n = 2 
(top) and n = 10 (bottom) in both pump modes. Dashed line 
shows Q parameter for generated fields, full line for pump 
fields. 
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FIG. 7: Mandel Q parameter for initial coherent states with 
(n) = 10 in both pump modes. Dashed line shows Q param- 
eter for generated fields, full line for pump fields. 



D. Realization of phase gate for 
continuous-variable quantum computation 

The field evolution from an initial Fock state with 
n = 1 photo ns in each of the pump modes, discussed 
in subsection [VB, exhibits an interesting feature. After 
a single complete cycle of energy conversion to the gen- 
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erated fields and back, the quantum state of the system 
undergoes a phase change of exactly 7r. This can be used 
to implement a phase gate for continuous quantum com- 
putation. Noting that no conversion occurs unless both 
input modes are excited and fixing the medium length at 
a value corresponding to twice the conversion length for 
a single photon input of both pump modes one has the 
following evolution of states 



10,0,0,0) 
|1, 0,0,0) 
10,1,0,0) 
|1.1,0,0) 



10,0,0,0), 
|1, 0,0,0), 
10,1,0,0), 
-|1, 1,0,0). 



(21) 



This is a perfect realization of a phase gate. 



V. MEAN-FIELD THEORY OF RESONANT 
FORWARD FOUR- WAVE MIXING 

In order to obtain a better understanding of the quan- 
tum dynamics of resonant four- wave mixing, in particular 
the limitations and scaling of conversion efficiency as well 
as the conversion distance, we now derive approximate 
analytical solutions within an appropriate mean-field the- 
ory. Our starting point is again the effective Hamiltonian 
(|l8|), along with the assumption that we can replace the 
denominator, a constant of motion, with its expectation 
value. In this case the equations of motion for the fields 
in the Heisenberg picture are 



di h = 


-— b 2 a ia2 


(22) 


d 

dt al = 


"Ad aA62 


(23) 




IKC „ „ 

"Ad haia2 


(24) 


d 

dT 2 = 


~Ad 016162 


(25) 



where d — (b\bi + djdi). 

The equations of motion for the average field intensi- 
ties (b[bi) and so on will contain four-field correlation 
functions, for example 



kt£t\ 



— {b\bi) = — {bib 2 a\d\) - (ciid 2 b\b : 



(26) 



where we have once again switched to a spatial rather 
than a temporal picture. 

To proceed further we assume that all fields are gaus- 
sian and that correlations between pump and generated 
fields can be neglected. This assumption is reasonable 
as long as the pump fields are initially in coherent states 
with sufficiently large amplitude. With this decorrelation 
approximation, the fourth order expectation values can 
be expressed as sets of bilinear terms of the form (b\bj) 
and (a\dj). It is important that we also keep anoma- 
lous correlations, such as (didj), as both pairs of fields 



are strongly correlated. In addition, we can use the con- 
stants of motion (^) - (||) to relate the expectation values 
of the number operators (b\b\), (blb 2 ) 7 (a\a\) and (a^az). 
If we define (6162) = bi2e lipb ', (0,10,2) = ai2e liPa , and 
(b\bi) = b, then we obtain the following set of coupled, 
complex- valued nonlinear differential equations: 

d 



dC ai2 



dc ai2 

He 6 



d^ a 



— (2a + l)ai2 sin(tp a - ip b ) (27) 

-£-(2b-2d-l)bi 2 sm(<p a -<p b ) (28) 
2k 



— b i2 E 12 sin(</j - ip b ) 



(29) 



-^(U+l^cos&a-w) (30) 
Ad 012 

-£-(2b-2d-l)-2-coB(<p a -<p b ). (31) 
Ad 012 



To simplify this system we note that if we take expecta- 
tion values of both sides of (0) we find b\ 2 a 12 cos((^ a — tp b ) 
is a constant of motion. If we assume that both gener- 
ated fields start from vacuum this constant is equal to 
zero, and since 612 and a\ 2 are not always zero we see 
that cos(</? Q — ipb) = ±1, with the sign flip coming at the 
end of each conversion cycle. This observation, in con- 
junction with two more constants of motion that can be 
extracted from (|^) - (|3l]), enable us to reduce the set of 
five coupled equations to just one: 



d& 



2 k 



^ = ^- Q V(b 2 +b-b Q )(b-b )(b-b -l), (32) 

where b = (SJ(O)Si(O)) = (SJ(0)S 2 (0)). This is the equa- 
tion governing the evolution of the expectation value of 
the number of photons in each of the pump modes, un- 
der the assumption that both pump fields have the same 
intensity and the two generated fields start from vacuum. 

As a fourth order polynomial is involved, this differ- 
ential equation can be (implicitly) solved analytically in 
terms of elliptic integrals, but the specific form of the so- 
lution is involved and not particularly illuminating. Nu- 
merical solutions to (|3^) are plotted in Figure || for a 
number of initial pump field intensitities. One recognizes 
an oscillatory exchange of energy with a non-perfect con- 
version effiency. Comparison of the mean-field results 
with the fully quantum calculation for the case of a co- 
herent pump input shows good agreement over the first 
period of energy exchange between pump and generated 
fields (Figure |9j). In particular the maximum conver- 
sion efficiency and the conversion distance are well re- 
produced. As in the quantum solutions, the conversion 
period increases logarithmically with the input intensity. 
For larger distances the mean-field solution remains pe- 
riodic, while the oscillations in the true quantum case 
decay. As the interaction distance increases, higher or- 
der correlations build up and thus the gaussian factoriza- 
tion approximation used in the mean-field theory breaks 
down. 
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FIG. 8: Pump field intensity in the mean-field approximation. 
From top to bottom: 10, 50, 100, 500 and 1000 photons. 
k/A = l,b = 10. 
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FIG. 9: Comparison between the exact numerical solutions 
(dashed line) and the mean-field theory (solid line) for a co- 
herent state with an average photon number of 1000. 



Analoguous to the scmiclassical case, ( B2[ ) can be 
mapped to a nonlinear pendulum problem, with b cor- 
responding to the position of a particle with mass A 2 /k 2 
moving in a potential 



v ( b ) = -72 ( &2 + b - M(& - b )(b -b -I) (33) 

where b plays the role of the particle position and the 
identification £ — > t is used. 



The potential (33) is plotted in Figure 10 for the case 
&o = 10, k/A = 1. From the shape of the potential 
one can see that different dynamics is expected com- 



/ 1000 








/ 500 








/ -500 




\ 5 




/ -1000 









FIG. 10: Nonlinear pendulum potential experienced by a 
"particle" whose position corresponds to the number of pho- 
tons in the pump modes. k/A = 1, bo = 10. 



pared with the semiclassical case, where the potential 
is described by Figure ||. In the semiclassical analysis, 
choosing the initial intensity of a\ and to be zero cor- 
responds to starting exactly on the right hand peak of 
the potential in Figure and consequently no dynami- 
cal evolution can occur. 

In the mean-field theory, however, the starting posi- 
tion (bo) is always slightly to the left of the right summit 
(peak position is ~ bo + | for large bo), and consequently 
oscillation will always take place. As the input power is 
increased the starting position asymptotically approaches 
the right hand peak, but never reaches it. From the po- 
tential ([53j) one can easily obtain the conversion distance 
or oscillation period: 



d& 



60 J- 2 -£v(b) 



(34) 



where 61 is the inner turning point. As the input power 
is increased, it takes longer for the oscillation to begin 
due to the flatter gradient of the potential, leading to 
the logarithmic increase in conversion period. This is in 
sharp contrast to ordinary off-resonance four-wave mix- 
ing, where the conversion length decreases with input 
power. The effective hamiltonian discussed in the present 
paper can be converted into that of ordinary off-resonant 
four-wave mixing by replacing the intensity dependent 
denominator d by a constant. Figure [ll] gives a com- 
parision of the conversion distance as a function of input 
power for the two cases. The periods were found from 
numerical solutions to the full quantum problem using 
the Hamiltonian (|3|) with and without the resonant de- 
nominator. One clearly recognizes the peculiar feature 
of the resonant process to work better for smaller input 
intensities. 

From the roots of (^) one can immediately obtain a re- 
sult connecting input intensity and conversion efficiency. 
Defining the conversion efficiency e = 6 m in/^o we find 



1 



e = l- 5 7-(v / 4V 



1-1) 



1 



1 - — . (35) 



Thus in the mean-field model the conversion efficiency 

— 1/2 

(1 — e) scales with the input pump field as b . 
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Input power 

FIG. 11: Scaling of conversion distance as function of input 
intensity of coherent pump for ordinary off-resonance four- 
wave mixing (dashed line, d — const.) and resonant four-wave 
mixing (full line). 



VI. SUMMARY 

In the present paper we have discussed the quantum 
theory of stationary resonant four-wave mixing in a for- 
ward scattering configuration. We considered the in- 
teraction of four electromagentic fields with atoms in a 
modified double-A configuration. The modified five-state 
coupling scheme was used instead of the standard four- 
level one in order to eliminate ac-Stark induced nonlinear 
phase shifts |27j]. This is necessary, because a full quan- 
tum analysis shows that full conversion cannot be ob- 
tained when these terms are present, even in the asymp- 
totic limit — 50% conversion is the best that can be 
achieved using coherent states. 

To eliminate the atomic variables we assumed adia- 
batic following and derived a classical nonlinear interac- 
tion Hamiltonian of the four electromagnetic fields. In 
contrast to ordinary off-resonant four-wave mixing, the 
denominator of this Hamiltonian contains the sum of the 
intensities of the resonant pair of pump and generated 
fields. This is a consequence of the infinitely long-lived 
two-photon resonance which is entirely determined by 
power-broadening for arbitrarily small field intensities. 
As a consequence the effective Hamiltonian cannot be ex- 
panded as a power series in the fields and the nonlinear 
interaction becomes stronger the smaller the pump field 
intensity. The classical evolution can be mapped onto a 
one-dimensional nonlinear pendulum and solved exactly, 
provided small seed fields are included in the analysis. 
The initial evolution of the generated fields from vacuum, 
however, is entirely determined by quantum fluctuations 
and thus cannot be described within a classical approach. 



To obtain a full quantum solution under stationary 
conditions we quantized the effective interaction Hamil- 
tonian by replacing c-number expressions by properly or- 
dered operator products. We showed that the constants 
of motion of the Hamiltonian allow us to reduce the sta- 
tionary wave mixing process to a single-mode problem, 
which can be solved numerically for up to 10 3 input pho- 
tons. We have analyzed the field evolution starting from 
initial Fock and coherent states in the pump modes. In 
both cases an oscillatory energy exchange between the 
modes was found. For the case of initial Fock states with 
up to n = 2 photons in both pump modes, the oscil- 
lation is purely sinusoidal with a single frequency and 
complete conversion is possible. For higher input pho- 
ton numbers multiple frequencies appear and complete 
conversion is only possible after several oscillation peri- 
ods. In the case of coherent input fields the oscillations 
in the energy exchange are damped and complete conver- 
sion can be achieved only asymptotically for large input 
intensities. 

As the four-wave mixing process simultaneously cre- 
ates photons in two modes, the photon numbers in those 
modes are perfectly quantum correlated, which has po- 
tential applications in quantum communication systems. 
The statistics of the photon pair production is, however, 
mostly super-Poissonian. 

We observed an interesting property of the four-wave 
mixing process when the two input fields are both in a 
single-photon Fock state. After a complete conversion 
cycle into the generated fields and back the quantum 
state attains a phase shift of exactly it. This allows one 
to construct an ideal phase gate for continuous variable 
quantum computation. 

In order to gain more physical insight in the quantum 
four-wave mixing we developed a mean-field theory as- 
suming gaussian and independent pump and generated 
fields but taking into account anomalous field correla- 
tions. The mean-field equations can again be mapped 
to a one-dimensional anharmonic pendulum and solved 
exactly. The solutions for the average intensities repro- 
duce the quantum results for coherent input over the first 
oscillation period. In particular analytic expressions for 
the conversion length and the conversion efficiency can be 
derived which coincide very well with the exact results. 
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